First-principles investigation of hyperfine interactions 
for nuclear spin entanglement in photo-excited fullerenes 



O 
(N 

1—5 

(N 



Vasileia Filidou,^'^ Davide Ceresoli,^'^ John J. L. Morton,^' ^ and Feliciano Giustino^ 

^Department of Materials, University of Oxford, Parks Road, Oxford 0X1 3PH, United Kingdom 
^CAESR, Clarendon Laboratory, Oxford University, Oxford, 0X1 3PU, United Kingdom 
^CNR-ISTM, via Golgi 19, 20133 Milan, Raly 
(Dated: January 30, 2012) 

The study of hyperfine interactions in optically excited fullerenes has recently acquired importance 
within the context of nuclear spin entanglement for quantum information technology. We here report 
a first-principles pseudopotential study of the hyperfine coupling parameters of optically excited 
fullerene derivatives as well as small organic radicals. The calculations are performed within the 
gauge- invariant projector-augmented wave method [C. Pickard and F. Mauri, Phys. Rev. B. 63, 
245101 (2001)]. In order to establish the accuracy of this methodology we compare our results 
with all-electron calculations and with experiment. In the case of fullerene derivatives we study 
the hyperfine coupling in the spin-triplet exciton state and compare our calculations with recent 
electron paramagnetic resonance measurements [M. Schaffry et aL, Phys. Rev. Lett. 104, 200501 
(2010)]. We discuss our results in light of a recent proposal for entangling remote nuclear spins in 
photo-excited chromophores. 
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I. INTRODUCTION 

In the area of quantum information processing sub- 
stantial efforts are being devoted to exploit the magnetic 
properties of molecular systems in order to realize quan- 
tum logic operations.^ In a recent proposal for entanghng 
two nuclear spins using a photo-excited transient electron 
spin^ the strength of the hyperfine interaction (HFI) has 
been identified as a key parameter determining the entan- 
glement efficiency. These investigations raise the ques- 
tion on how to optimize such parameters, and how to 
design new materials with tailored magnetic properties 
for specific quantum logic operations. 

First-principles calculations based on density- 
functional theory (DFT) have recently met with a 
number of successes in the calculation of the magnetic 
properties of materials, in particular when it comes to 
the building blocks of the spin Hamiltonian such as 
chemical shifts^"— J- couplings ^2J£ and g- tensor sJ^ii^ 
As magnetic spin resonance experimental techniques 
become more accurate, and novel applications to physics, 
biology, and materials science are being explored^ii^ 
DFT-based modelling of charge/spin interactions hold 
the promise of becoming a key player in the interpreta- 
tion of measured magnetic responses and in the design 
of materials with tailored functionality. 

In this work we aim to establish the validity and accu- 
racy of DFT-based modelling in the calculation of HFI 
parameters in a class of materials of interest for quantum 
information technology. We here focus on the HFI inter- 
actions in the fullerene derivatives synthesized in Ref. [H, 
since electron spin resonance measurements on these ma- 
terials have already been performed. 

The manuscript is organized as follows. In Sec. HIl we 
briefly review the connection between the HFI parame- 
ters and the efficiency of the quantum entanglement op- 



erations introduced in Ref. '2". In Sec. [1111 we describe 
the computational methodology that we employ for cal- 
culating hyperfine tensors. In Sec. IIV Al we discuss our 
HFI calculations for a series of small radicals. In this 
section we establish the accuracy and reliability of pseu- 
dopotential (PP) calculations of the HFI by comparing 
with all-electron (AE) calculations and with experiment. 
In Sec. IIVBI we present our calculations of the HFI on 
photo-excited fullerene derivatives, and we discuss our 
results in light of the entanglement protocol proposed in 
Ref. 0. In Sec. El we present our conclusions. 



II. ENTANGLING NUCLEAR SPINS VIA THE 
HYPERFINE INTERACTION 

In Ref. [2I a quantum information processing protocol 
involving nuclear spins has been introduced. This proto- 
col involves the interaction of two nuclear spins attached 
to a fullerene molecule through functional additions car- 
rying isotopically labeled nuclei. The ground state of 
this fullerene derivative is spin-unpolarized, and hence 
the HFI of the nuclear spins is vanishing. Upon optical 
excitation a spin-triplet excitonic state is generated in the 
molecule, and the nuclear spins interact with the spin- 
polarized electrons leading to a non- vanishing hyperfine 
coupling. As a result the two nuclear spins experience 
an interaction which is mediated by the spin of the exci- 
ton. Such interaction can be used to perform an entan- 
glement operation between the two nuclear spins. When 
the exciton decays and the fullerene reverts to its spin- 
unpolarized ground state, the nuclear spins may remain 
entangled. The efficiency of this protocol results from 
the competition between three timescales: (i) the time 
required to perform the entanglement operation, (ii) the 
exciton lifetime, and (iii) the decoherence time of the en- 
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tangled state The competition between these timescales 
defines an optimum range of HFI strength for achieving 
maximum entanglement. 

In order to establish the feasibility of the protocol pro- 
posed in Ref. 2, pulsed electron paramagnetic resonance 
experiments were performed on a prototype system, the 
fullerene diethyl malonate molecule [C60-C7H12O4, sys- 
tem Ml in Fig. [U^a)]. Although this system carries only 
one isotopically labeled nuclear spin associated with the 
^^C atom of the functional group directly bonded to the 
cage, it was possible to extract relevant information such 
as the lifetimes of the exciton states and the HFI param- 
eters. For the system considered the measured relaxation 
rates of 90-100 /is at 20-50 K correspond to maximum en- 
tangling power when the HFI parameters are in the range 
2 MHz to 12 MHz, with an optimal value of 6 MHz.^ The 
measured HFI parameter of 3 MHz in this prototype sys- 
tem with a single nuclear spin falls within the optimum 
theoretical range and motivates further work in this area. 

The key open question which remains is how to identify 
systems involving multiple nuclear spins and exhibiting 
optimal HFI strengths. One possible route for answer- 
ing such question is to employ first-principles materials 
modelling based on density-functional theory in order to 
screen candidate systems prior to their experimental syn- 
thesis. Accordingly, in the present work we aim at estab- 
lishing the accuracy of such first-principles calculations, 
in view of optimizing and designing functional materials 
for spin-based quantum information processing. 



to the nuclear site and r = |r|. The component a^^^ 
of the tensor A^^^ provides the isotropic hyperfine in- 
teraction and is referred to as the Fermi contact term. 
The traceless tensor h^^^ gives the anisotropic or dipo- 
lar hyperfine interaction. The anisotropic tensor h^^^ is 
conventionally given in terms of its eigenvaules 6^, ^yy^ 
and }}zz ' We note that the principal axis of the hyperfine 
tensor A^^^ and of its anisotropic part h^^^ do coincide by 
construction. 

The evaluation of the HFI parameters using Eqs. ([2]) 
and (|3]) requires the calculation of the electron spin den- 
sity. In the dipole-dipole term h^^^ the electron spin- 
density is weighted by the factor and the contribu- 
tions outside of the core region are dominant. In the 
Fermi contact term the electron spin density must be 
evaluated at the position of the nucleus, therefore an 
accurate description of the electron density within the 
core region is essential. While standard pseudopotential 
methods can be used for calculating the dipole-dipole 
term h^^\ the calculation of the Fermi contact contri- 
bution requires the use of all-electron methods or the 
projected augmented wave (PAW) method.— The PAW 
method enables the reconstruction of the all-electron 
wavefunction inside the core region, starting from a pseu- 
dopotential calculation of the valence electronic struc- 
ture, and has successfully been applied to the calculation 
of hyperfine parameters 



III. COMPUTATIONAL METHODOLOGY 

A. Hyperfine tensor from first principles 

The coupling of the electron spin 5 to a set of N nu- 
clear spins Ii {I = 1^...N) can be described using the 
following hyperfine Hamiltonian)^ 

H= 5^AW/1, (1) 

1=1, N 

A*^^^ being the hyperfine coupling tensor associated with 
the l-th nucleus at the site In atomic units the hy- 
perfine tensor can be written as)^ A^^j = a^^^dij + 
with 

tt^^^ = -j9el^egil^Nps{^l)^ (2) 

and 

% =gefJ^egH^N / dr — ps(r). (3) 

In these equations Ps(r) = Pt(^) ~ Pi(^) electron 
spin density, the free-electron g-factor, /ig the Bohr 
magneton, gi the gyromagnetic ratio of the l-th nucleus 
and fiN the nuclear magneton. In Eq. (|3]) r is relative 



B. Systems considered 

In Ref. Q we reported a preliminary calculation of the 
isotropic component of the HFI tensor of the Ml system 
[Fig. [T]^a)], and of a related fullerene derivative carrying 
two ^^C nuclear spins [system Bl in Fig. [TJb)]. Here we 
complement the work of Ref. ^ by calculating the entire 
HFI tensor of these systems, and we further extend it 
to fullerene derivatives of various conformations. Fur- 
thermore, we perform calculations on a number of small 
open-shell molecules, both within PP and AE methods, 
in order to establish the accuracy of the computational 
methodology. 

We consider following small molecules containing hy- 
drogen and oxygen: CH2CH3, CH2CH, CO3, H3CO, 
HCO, H2CO+, and CH3. In the case of diethyl mal- 
onate fullerene derivatives we consider the mono-adduct 
Ml [Fig. [B^a)], as well as the three regioisomeric bis- 
adducts B1-B3 shown in Fig. [T]^b). During the synthesis 
of stable regioisomeric bis-adducts from diethyl malonate 
and fullerene molecules, the trans- 1 isomer [model Bl in 
Fig- [11(b)] is produced with the lowest yield (0.8%), while 
the equatorial and the trans-3 isomers [models B2 and B3 
in Fig. [TJb), respectively] are produced with the highest 
yields (16% and 12%, respectively ) 
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B1 B2 B3 

(trans-1) (equatorial) (trans-3) 



FIG. 1: (a) Ball-and-stick (left) and schematic (right) rep- 
resentations of the diethyl malonate (C7H12O4) Ceo mono- 
adduct (Ml). The bridging ^^C atom carrying a nuclear spin 
is indicated. In the schematic representation the red bond 
indicates the anchoring point of the adduct on the fullerene 
cage. Atomic color code: C (grey), O (red), H (white), (b) 
Schematic representations of three regioisomeric diethyl mal- 
onate fullerene bis-adducts (B1-B3). 



C. Calculations details 

We perform pseudopotential (PP) density-functional 
theory calculations within the generalized gradient ap- 
proximation of Perdew, Burke, and Ernzerhof (PBE).^^ 
We describe the valence electronic wavefunctions using a 
planewaves basis set with a kinetic energy cutoff of 70 Ry. 
We take into account the core-valence interaction using 
Troulier-Martins normconserving pseudopotentials^i^i^S 
For the calculation of the HFI tensor we use the gauge- 
including projector- augmented wave method (GIPAW) 
of Ref. [2l|. We use periodic simulation cells, and in or- 
der to minimize interactions between periodic replicas of 

o 

the molecules, we use cells of size 12 A for the small 

o 

molecules, and cells of size up to 24 A for the fullerene 
derivatives. All PP calculations are performed using the 
quantum-ESPRESSO software package. 

For the calculation of the Fermi contact term, the ac- 
curate description of the electronic wavefunctions in the 
core region is required. In pseudopotential methods only 
valence electrons are described, and in some cases this 
might lead to significant errors due to the neglect of core 
polarization effects. Core polarization refers to the spin 
polarization of core electrons induced by the spin of va- 
lence electrons, and leads to an additional contribution 
to the Fermi contact term j^^i^^ Within PP calculations 



we introduce core polarization corrections using the per- 
turbative approach of Ref. ^ as follows. For each atom 
with core electrons (C and O) we calculate the spher- 
ical average of the spin density around the nuclues up 
to a radius of 5 a.u. using a logarithmic radial grid. 
This average is used to evaluate the perturbation AV 
due to valence polarization following Eq. (20) of Ref. [23 , 
and subsequently the contribution to the Fermi contact 
term arising from the core polarization using Eq. (17) of 
Ref. [23. The atomic all-electron wavefunctions are stored 
in the pseudopotentials files and the integrations are car- 
ried out on the radial grid. In the following, pseudopo- 
tential calculations including the core polarization cor- 
rection will be referred as "PP-CP" , while those without 
core polarization will be denoted simply by "PP" . In the 
case of light atoms such as those considered in this work 
relativistic effects are expected to be negligible, therefore 
we generated the pseudopotentials by solving the non- 
relativistic Schrodinger equation (it is custom practice 
to solve the scalar-relativistic equation only for atoms 
starting from the fourth row of the periodic table). In 
evaluating the Fermi contact term, the first relativistic 
correction can be taken into account through a spherical 
integration of the spin density up to the Thomson ra- 
dius tt, which provides an estimate of the finite size of 
the nucleus. However, the Thomson radii of H and C fall 
around the first point of our pseudopotential radial grid, 
therefore it seems more appropriate to evaluate the spin 
density at the nucleus (this is done by extrapolating the 
three innermost values on the grid). 

In order to establish the accuracy of PP calculations of 
HFI parameters we also perform AE calculations using 
the PBE exchange and correlation, as implemented in the 
Gaussian 09 program package For the AE calculation 
of the HFI tensor we adopt the EPR-II basis, which has 
been purposely designed to accurately describe the Fermi 
contact term^^i 

In order to study the HFI in optically-excited fullerene 
derivatives we perform calculations using the ground 
state of the Sz = I electron spin configuration. In fact, 
while density-functional theory does not describe cor- 
rectly excited states, the lowest spin-triplet state of Ceo, 
with the HOMO and LUMO states singly occupied with 
parallel spins, is the ground state of the Sz = I spin 
configuration, hence the use of DFT is legitimate 

For each structure M1-B3 shown in Fig.[T]we optimize 
several initial geometries in both the singlet Sz = and 
triplet Sz = I configurations in order to identify the cor- 
responding ground state. The energy difference between 
the triplet and the singlet states ranges from 1.29 to 1.33 
eV across all the models considered. This result is in 
agreement with the measured excitation energy of 1.30 
eV of the non-functionalized fullerene Ceo'^'^ This sug- 
gests that our description of the lowest triplet exciton 
state is adequate. 
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TABLE I: Fermi contact hyperfine coupling parameters of 
small molecules containing C, and H. We report calcula- 
tions using the GIPAW method without accounting for core- 
polarization effects (PP), calculations within GIPAW includ- 
ing core-polarization corrections (PP-CP), all-electron calcu- 
lations (AE), and measured Fermi contact parameters. The 
atoms are listed according to corresponding chemical formula. 



Molecule Atom Fermi contact term (MHz) 







PP 


PP-CP 


AE 


Expt. 


CH2CH 


C(l) 


397 


305 


306 


302^ 




H(l) 


125 




124 


96 




H(2) 


191 




190 


192 




C(2) 


-27 


-17 


-18 


-24 




H(3) 


38 




43 


37 


CH2CH3 


C 


205 


81 


79 


83^ 




H 


-63 




-61 


-69 




C 


-40 


-38 


-34 


-38 




H 


51 




40 


75 


CHs 


C 


214 


79 


72 


80 s 




H 


-69 




-66 


-69 


H3CO 


H 


94 




95 


122 




C 


-35 


-44 


-37 


_44d 


HCO 


H 


379 




367 


356 




C 


435 


428 


389 


367^ 


H2CO+ 


H 


279 




372 


372 




C 


-76 


-83 


-87 


-109^ 


CO3 


C 


-31 


-31 


-28 


-32^ 



^Refsj3p|l3l|, 
^Ref. 
"Ref. I 
^Ref. 
^Ref. 
^Ref. 
SRef. 



IV. RESULTS AND DISCUSSION 
A. HFI parameters of small molecules 

In Table |I] we present the results of our pseudopoten- 
tial calculations and we compare with all-electron calcu- 
lations and with experiment. For a quantitative analysis 
we present both pseudopotential results with (PP-CP) 
and without (PP) core polarization correction. 

In the cases of CH2CH, CH2CH3, and CH3 Table |I] 
shows that the neglect of core polarization for the C 
atoms leads to an underestimation of the Fermi contact 
term by as much as 92-135 MHz. When the core polar- 
ization correction is introduced, all the calculations tend 
to agree rather well with our AE results and with experi- 
ment. In the case of H atoms, where the core and valence 
do coincide, our PP results are in good agreement with 
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FIG. 2: Comparison between the HFI parameters calculated 
using PP-CP and AE methods and experimental data, from 
Tables HI and [III (a) Fermi contact term of C atoms, (b) Fermi 
contact term of H atoms, (c) Dipolar coupling of C atoms, 
(d) Dipolar coupling of H atoms. If the calculations were able 
to predict experimental data exactly, then all the data points 
would fall on the diagonal solid line. 



AE and experimental data. 

Figures [2]^a),(b) show a comparison between the Fermi 
contact term for C and H atoms calculated with the PP- 
CP method or the AE method, and experiment. The 
relative accuracy of PP-CP calculations and AE calcu- 
lations in the case of C atoms appears similar for this 
set of test molecules. Indeed the largest deviations from 
experiment are of 17% and 13% for PP-CP and AE cal- 
culations, respectively. In absolute terms the largest de- 
viations from experiment are 61 MHz (22 MHz) for C 
atoms using PP-CP (AE) calculations, and 93 MHz (28 
MHz) for H atoms using PP-CP (AE) calculations. 

In Table [III we present the results of calculated dipo- 
lar hyperfine couplings. The same data set is reported 
in Figs. [2fc),(d) for C and H atoms, respectively. The 
largest absolute deviation between our PP calculations 
and experiment is 23 MHz, while the largest deviation 
for the case of AE calculations is 21 MHz. Again PP 
calculations and AE calculations appear to perform sim- 
ilarly, and are both able to provide reasonably reliable 
data when compared to experiment. 

What we learn from these tests is that, on the one 
hand, in the case of the Fermi contact term, the inclu- 
sion of core polarization in PP calculations is critical to 
achieving agreement with AE calculations and with ex- 
periment. On the other hand, PP calculations of the 
dipolar HFI interaction are already in reasonable agree- 
ment with AE data and experiment. In the case of HCO 
we observe that the deviation between the calculated 
isotropic HFI (either PP-CP or AE) and experiment can 
be rather large, hence care must be used in the interpreta- 
tion of the calculated HFI parameters. This discrepancy 
requires further investigation. 
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TABLE II: Dipolar hyperfine coupling parameters of small molecules containing C and H. We report PP calculations, AE 
calculations, and measured parameters. The atoms are listed according to corresponding chemical formula. 



Molecule Atom Dipolar HFI (MHz) 


Molecule Atom Dipolar HFI (MHz) 


Axis PP AE Exp. 


Axis PW AE Exp. 



CH2CH 




bxx 


-86 


-72 


H3CO 




bxx 


7 


6 






C(l) 


byy 
b.. 


-64 
148 


-56 
127 




C 


byy 

bzz 


2 
-9 


2 
-8 










37 


37 


38^ 




bxx 


-8 


-7 






H(3) 


byy 


-6 


-6 


-11 


H 


byy 


-4 


-4 










-32 


-32 


-27 




&ZZ 


12 


12 




CH2CH3 




bxx 


-83 


-73 


HCO 




bxx 


-43 


-42 


-48^ 




C 


byy 


-82 


-72 




C 


byy 


-33 


-38 


-24 






b,. 


165 


145 






bzz 


76 


80 


72 






bxx 


35 


38 


35^ 




bxx 


-14 


-14 


-14 




H 


byy 


2 





1 


H 


byy 


-8 


-7 


-8 






b,. 


-38 


-38 


-36 




bzz 


23 


21 


23 






bxx 


-3 


-3 


H2CO+ 




bxx 


-25 


-21 


-15^ 




C 


byy 


1 


1 




C 


byy 


7 


6 


4 






6.. 


3 


2 






6zz 


18 


16 


10 






bxx 


-4 


-4 


-3^ 




bxx 


-10 


-10 


-9 




H 


byy 


-3 


-3 


-3 


H 


byy 


-8 


-7 


5 






b.. 


7 


9 


6 




6zz 


18 


18 


4 


CH3 




bxx 


-70 


-74 


-63^ CO3 




bxx 


5 


4 


4j 




C 


byy 


-72 


-74 


-63 


C 


byy 


5 


4 


4 






6.. 


150 


148 


127 




bzz 


-9 


-9 


-8 






bxx 


-36 


-39 


-36"^ 














H 


byy 

bzz 


-1 
35 


-1 
39 


-1 

35 













Ref . S ' Ref . 113, ' Ref . S ^ Ref . [H, ^ Ref . Ell, Ref . E, " Ref . S. 



B. HFI parameters of fuUerene derivatives 

In this section we present the HFI parameters calcu- 
lated for the fullerene derivatives of Fig. [H Figure [3] 
shows the calculated spin density distribution for the 
structures Ml and Bl-3 in the triplet state. The spin den- 
sity is greatest (blue regions) on the fullerene cages, while 
it is rather small (red regions) on the adducts and the 
isotopically labeled C nuclei. As the strength of the HFI 
scales with the spin density [cf. Eqs. (O,©], Fig. [3] sug- 
gests that we should expect a rather small HFI strength 
for all the structures considered. 

This observation is consistent with the calculated HFI 
parameters reported in Table IIIH which are of the or- 
der of only a few MHz. All the pseudopotential calcula- 
tions in this table include the core polarization correction 
(PP-CP) discussed in Sec. HITCl For the case of the Ml 
structure both PP-CP and AE calculations are in_good 
agreement with the experimental data from Ref. [2|. In 
light of the error bars established in Sec. IIV Al for the 
test molecules, this level of agreement (around 1 MHz) 




FIG. 3: Spin density distribution in the lowest triplet excited 
state of the four fullerene derivatives Ml and Bl-3 considered 
in this work. The color code indicates the magnitude of the 
spin density (the density isosurface corresponds to a isovalue 

of 0.07 electron/A^). 
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top view side view 



FIG. 4: Ball-and-stick representation of the diethyl malonate 
adduct on a portion of the Ceo fullerene (color code as in 
Fig. [T]). The angles in all the systems considered Ml and 
Bl-3 differ by less than one degree: a, a =126°, /3, /3'=245°, 
7, 7^= 111°, ^=112°, 6>=65°. The different C-C bond length 
dl in the two adducts of B2 lifts the symmetry between the 
two isotopically labeled nuclei and results into different HFI 
parameters [cf. B2(l) and B2(2) in Table Im]. 



may be somewhat fortuitous. However it is reasonable 
to assume that, due to error cancellation, our calcula- 
tions should provide a reliable description of the relative 
magnitude of the HFI parameters across various fullerene 
derivatives. It is possible that accounting for thermal mo- 
tion may lead to slightly different HFI parameters in our 
calculations, however we expect that such effects will fall 
within our error bars established in Sec. IIV A[ Morover 
the effect of thermal motion should be similar across the 
various derivatives, therefore the relative magnitude of 
the HFI parameters is expected to be reliable. 

It is interesting to observe in Table Hill that the location 
of the diethyl malonate adduct has a significant effect on 
the HFI parameters, leading to variations of up to an 
order of magnitude in the Fermi contact term. 

For both the trans- 1 isomer Bl and the trans-3 isomer 
B3 Table mil shows that the two ^^C nuclei have the same 
HFI tensor. At variance with this, in the equatorial con- 
former B2 only one nucleus exhibits a sizeable HFI. In 
order to rationalize this finding we examine the structural 
parameters of each fullerene derivative. Figure. 21 shows 
the geometry of the diethyl malonate adducts on the Ceo 
cage. The structure of the adducts (i.e. bond lengths and 
bond angles) are almost identical in each fullerene deriva- 
tive. The only noticeable difference is that the length 
of the C-C bond bridging the adduct to the fullerene is 

o 

1.58 A for both adducts in the structures Bl and B3, 
while in the structure B2 the two adducts have different 

o 

C-C bond lengths of 1.60 and 1.57 A, respectively. This 
indicates that in the equatorial conformer B2 one of the 
two nuclear spins is closer than the other to the Ceo cage. 
Since the spin density is largest on the fullerene cage, the 
nuclear spin with the shorter C-C bond length can be ex- 
pected to exhibit a larger HFI strength, consistent with 
the results in Table Hill 

In Ref. i it was shown that, in systems like the struc- 
tures Bl-3 considered here, it should be possible to 



TABLE III: Calculated HFI parameters of the fullerene 
derivatives Ml and Bl-3 shown in Fig. [1] We present data 
from pseudopotential calculations including core polarization 
correction (PP-CP), and we also provide the calculated Fermi 
contact term without core polarization correction (PP) for 
completeness. The dipolar terms are the same in both cases. 
For the fullerene mono- adduct Ml we also provide AE data 
and experimental data from Ref. IH (measured at 20 K) for 
comparison. In the case of the Bl-3 systems the values for 
the two ^^C-labeled nuclei are reported as (1) and (2), respec- 
tively. All the data are in MHz. 



Structure 


a 


bccx 


^yy 


b.. 


Ml, PP-CP 


2.0 


-1.6 


-2.3 


4.0 


Ml, PP 


4.1 








Ml, AE 


1.0 


-1.3 


-2.0 


3.4 


Ml, Expt. 


2.0 








Bl(l), PP-CP 


1.0 


-0.7 


-1.2 


1.9 


Bl(l), PP 


2.2 








Bl(2), PP-CP 


1.0 


-0.7 


-1.2 


1.9 


Bl(2), PP 


2.2 








B2(l), PP-CP 


0.2 


0.0 


0.6 


-0.6 


B2(l), PP 


0.3 








B2(2), PP-CP 


1.5 


-2.3 


-2.4 


4.7 


B2(2), PP 


4.4 








B3(l), PP-CP 


2.2 


-1.2 


-1.8 


3.1 


B3(l), PP 


3.9 








B3(2), PP-CP 


2.2 


-1.2 


-1.8 


3.1 


B3(2), PP 


3.9 









achieve the entanglement of the two nuclear spins if the 
following conditions are fulfilled: (i) the Fermi contact in- 
teraction when the system is photo-excited to its triplet 
state is in the range 2-16 MHz. (ii) the two isotopically 
labeled C nuclei exhibit a very similar strength of the 
HFI. Table HIH shows that, on the one hand, the trans-1 
Bl and the equatorial B2 isomers do not fulfill at least 
one of these requirements. On the other hand, the trans-3 
B3 isomer fulfills both criteria, and hence it represents a 
promising candidate for nuclear spin entanglement. This 
finding is particularly interesting since the B3 structure 
can be produced with the relatively high yield of 12% (cf. 
Sec HITBl) . 



V. CONCLUSIONS 

In this work we perform calculations of the hyper- 
fine couling tensor for photo-excited fullerene derivatives 
of interest in spin-based quantum information process- 
ing. We report extensive validation of pseudopotential 
methods against all-electron calculations and experimen- 
tal data for a series of small radicals. A key observation is 
that PP calculations can deliver an accuracy comparable 
to AE calculations, provided the effect of core polariza- 
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tion is included. 

We investigate the HFI strength in a series of photo- 
excited diethyl malonate fullerene derivates, which are 
potential candidates for nuclear spin entanglement oper- 
ations. We clarify the relation between the calculated 
HFI parameters and the underlying atomic-scale struc- 
ture of these systems. We identify one particular struc- 
ture, the trans-3 isomer B3, as the best match to the en- 
tanglement criteria set in Ref. IH. The fullerene derivative 
thus identified can be produced with relatively high yield 
and thus represents a potential candidate for achieving 
the entanglement of two nuclear spins in a photo-excited 
chromophore. 

Our work provides a set of validation tests and bench- 
marks for future calculations of HFI parameters using 
pseudopotential methods. In addition, this work repre- 
sents a first step towards the synergistic use of atomistic 
modelling and electron spin resonance experiments in the 



design of functional materials for quantum information 
technology. 
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